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ABSTRACT 

We present a systematic analysis of the intracluster medium (ICM) in an X-ray flux limited sample 
of 45 galaxy clusters. Using archival ROSAT PSPC data and published ICM temperatures, we present 
best fit double and single /3 model profiles, and extract ICM central densities and radial distributions. 
We use the data and an ensemble of numerical cluster simulations to quantify sources of uncertainty for 
all reported parameters. 

We examine the ensemble properties within the context of models of structure formation and feed- 
back from galactic winds. We present best fit ICM mass-temperature Mjcm-{Tx) relations for Mjcm 
calculated within r5oo and Ih^Q Mpc. These relations exhibit small scatter (17%), providing evidence 
of regularity in large, X-ray fiux limited cluster ensembles. Interestingly, the slope of the Mjcm-{Tx) 
relation (at limiting radius r^oo) is steeper than the self-similar expectation by 4.3cr. We show that there 
is a mild dependence of ICM mass fraction ficM on (Tx); the clusters with ICM temperatures below 
5 keV have a mean ICM mass fraction {ficAi) = 0.160 ± 0.008 which is significantly lower than that of 
the hotter clusters {ficM) = 0.212 ± 0.006 (90% confidence intervals). In apparent contradiction with 
previously published analyses, our large. X-ray fiux limited cluster sample provides no evidence for a 
more extended radial ICM distribution in low {Tx) clusters down to the sample limit of 2.4 keV. 

By analysing simulated clusters we find that density variations enhance the cluster X-ray emission 
and cause Mjcm and ficM to be overestimated by ~12%. Additionally, we use the simulations to 
estimate an ficM depletion factor at rsoo. We use the bias corrected mean ficM within the hotter 
cluster subsample as a lower limit on the cluster baryon fraction. In combination with nucleosynthesis 
constraints this measure provides a firm upper limit on the cosmological density parameter for clustered 
matter < (0.36 ± 0.01)/i5o^/^ 

Subject headings: galaxies: clusters: general — intergalactic medium — cosmology 



1. INTRODUCTION 

The properties of galaxy cluster virial regions provide 
powerful constraints on models of structure formation and 
evolution. For example, evidence for continuing accretion 
in nearby clusters constrains the cosmological density pa- 
rameter Oq and the slope of the power spectrum of den- 
sity perturbations P{k) on cluster scales. Combining the 
baryon fraction within cluster virial regions with nucle- 
osynthesis constraints on the baryon to photon ratio pro- 
vides an estimate of the cosmological density parameter 
for clustered matter flM- An emerging theoretical con- 
sensus that dissipationless collapse produ ces haloes which 
conform to a "universal" density profile ( Navarro, Frenk] 
& White 1997) implies that measurements of cluster virial 
structure provide constraints on the nature of dark mat- 
ter. Therefore, systematic, high precision analyses of large, 
well defined cluster samples can potentially enjoy broad 
and lasting impact. 

There are several tools available to probe galaxy cluster 
virial regions. One can use gravitational lensing to study 
the characteristics of the clustered mass directly. At the 
present time, typical cluster weak lensing maps have peak 
signal to noise of ~10 even at disappointingly low angular 
resolutions (e.g. Fischer & Mohr in prep), making it dif- 
ficult to resolve differences between observed profiles and 



the NEW prediction. Studies of strong lensing can provide 
extre mely detailed constraints on the structure of cluster 
cores (Tyson, Kochanski & DellAntonio 1998), but the 
required alignments of lens, source and observer occur in 
few clusters. Dynamical studies using galaxies can poten- 
tially resolve the dark matter structure and galactic orbits 
simultaneously, but doing so requires such a large galaxy 
sample that data from multiple clust ers must be rescaled 
and stacked (e.g. [Carlberg et al 19*97 ), precluding tests of 
anything except average properties of cluster mass profiles. 
Archival ROSAT PSPC X-ray images of nearby galaxy 
clusters are of sufficient quality that high precision studies 
of the ICM in individual clusters are possible. Moreover, 
it is straightforward to study a large, well defined clus- 
ter ensemble, making it possible to compare the structure 
of high and low mass clusters and to quantify departures 
from self similarity. 

Here we present a precision analysis of the intraclus- 
ter medium in 45 nearby clusters. Our goals include (i) a 
study of cluster regularity, (ii) a quantitative characteriza- 
tion of the ICM mass-temperature {Micm-{Tx)) relation 
for comparison to future structure formation and galaxy 
formation simulations, (iii) a definitive test of whether the 
ICM mass fraction ficM and radial distribution vary with 
cluster mass, and (iv) an improved measure of {fichi), 
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leading to an upper limit on the cosmological density pa- 
rameter rijv/- 

Our technique is a synthes is of profile fitting methods 



appli ed to ell iptical galaxies (Saglia et al. f993, Mohr & 



Wegn^r 1997) and procedu res applied to examine the ICM 
in sm aller cluste r samples ( David et al. 1993| , David, Jone^ 
& For nan 1995). The primary observational data required 
to study each cluster are (i) the X-ray surface brightness 
profile, which constrains the radial distribution of the ICM 
along with (ii) the X-ray luminosity within some aperture 
and (iii) the emission weighted mean ICM temperature 
(Tx), which together constrain the ICM central density. 
We demonstrate that our ICM constraints are insensitive 
to temperature variations within the cluster. Departing 
from some other published analyses, we do not assume a 
particular form for the dark matter distribution, because it 
is not necessary to do so when studying the ICM; however, 
we do assume a spherically symmetric ICM distribution. 
We test our technique on an ensemble of hydrodynami- 
cal cluster simulations and quantify the effects of present 
epoch mergers and asphericity on our results. 

We describe the data and numerical simulations in §2, 
detail the analysis in §3, and then present the results in 
§4. Section 5 contains a discussion of the results. We use 
Ho = 50 km/s/Mpc and Qo — ^ throughout. 

2. DATA 

We study members of an X-ray flux limited sample of 
55 clusters (Edge et al. 1990) that were observed with the 
ROSAT PSPC. Archival PSPC images of forty-eight mem- 
bers exist, but Abell399 lies partially outside the PSPC 
field of view (FOV), Virgo extends well beyond the PSPC 
FOV, and the quality of the short (850 s) exposure of 
Abell2147 image is well below the others. Of the remain- 
ing 45 clusters, three have no upper limits on their tem- 
perature uncertainties, and so as noted later we exclude 
them from portions of the final analysis. The final sam- 
ple is very similar to the same sample an alyzed in a study 



of th e X-ray Size- Temperature Relation ( Mohr 
'19971). 



Evrard 



2.1. X-ray Images 

We obtain the ROSAT data products on-line through 
the High Energy Astrophysics Science Archive Research 
Center (HEASARC), and we reduce the se images using 



PROS and Snowden analysis software (Snowden et al 
1994). The X-ray image for each PSPC observation is the 



sum of individually flat fielded Snowden energy bands i?4 
through i?7, corresponding approximately to the energy 
band 0.5-2.0 keV. We exclude time intervals with master 
veto rates higher than 220cts/s, and exclude other high 
background time intervals whose inclusion would degrade 
the detection significance of a source 10% as br ight as the 
background (Pildis, Bregman & Evrard 1995). We pro- 
duce uncertainty images by assuming Poisson noise, es- 
timating a 1% uncertainty in the flatfields, and applying 
standard error propagation. Multiple PSPC observations 
of the same cluster are registered and combined using the 
positions of bright point sources where possible, or, al- 
ternatively the image header pointing positions. We re- 
move obvious point sources not coincident with the peak 
in the extended emission using the IRAF task imedit. The 



pixel scale is 14'.' 947, and the effective angular resolution is 
FWIIM~0.'75. Table contains a Hst of clusters and rel- 
evant information, including ROSAT sequence number(s) 
of observations incorporated into the final cluster image, 
effective exposure time t^xp, and the emission weighted 
mean ICM temperature {Tx)- 

The procedure described above produces images in units 
of photons per second; we use PROS software to convert 
to cgs flux units ergs/s/cm^ in the restframe of the cluster. 
In the conversion we assume an isothermal ICM emitting 
a Raymond-Smith spectrum with 30% solar abundances 
and a temperature (Tx). We also correct for galactic ab- 
sorption. 
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Fig. 1.— We plot the relative PSPcfcount rate in the 0.5-2.0 keV 
band for a parcel of gas with constant emission measure and varying 
temperature. The emission model is Raymond-Smith v^ith 30% so- 
lar abundances. Note that for Tx > 1.5 keV a factor of two increase 
in temperature only decreases the PSPC count rate by ~10%. 

2.2. ICM Temperatures 

We use published, emission weighted ICM temperatures 
(Tx) in our analysis (see Table 0). Approximately a third 
of these temperatures come from Einstein MPC observa- 
tions, but where possible we have substituted more ac- 
curate measurements made with Ginga or ASCA obser- 
vations. The sources fo r the temperature measur ements 
include: Einstein MPC: [David et al. 1993: Ginga: Day el 



al. 1991, Allen et al. 1992, lohnstone et al. 1992, Davie 



et al. 19931 , jHughes et al. 1993|, [Arnaud fc Evrarc 



19981 ASCA: [Henriksen fc Markevitch 1996|, iMarkevitch 



1996, Matsuzawa et al. 1996, Tamura et al 199£, Marke 



vitch & Vikhlinin 1997, Markevitch et al. 1998 



2.3. Hydrodynamical Cluster Simulations 

We use an ensemble of 48 hydrodynamical cluster sim- 
ulations to better understand the systematics in the mea- 
surements presented below. The cluster simulations are 
carried out within four different cold dark matter (CDM) 
dominated cosmogonies (1) SCDM {fl = 1, erg = 0.6, 
h^o = 1, r = 0.5), (2) TCDM(r! = 1, as = 0.6, h^o = 1, 
r = 0.24), (3) OCDM {no = 0.3, as = 1.0, /iso = 1-6, 
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r = 0.24), and (4) LCDM (f^o = 0.3, Xq = 0.7, as = 1.0, 
h^o — 1.6, r = 0.24). Here CTg is the power spectrum 
normalization on 8h^^ Mpc scales; initial conditions are 
Gaussian random fields consi stent with a CDM transfer 
function with the specified T ( Davis et al. 1985| ). Within 
each of these models, we use two 128'^ N-body only simula- 
tions of cubic regions with scale ~ 400 Mpc to determine 
sites of cluster formation. Within these initial runs the 
virial regions of clusters with Coma-like masses of 10^^ Mq 
contain ~10"^ particles. 

Using the N-body results for each model, we choose 
clusters for additional study. We zoom in on these clusters, 
resimulating them at higher resolution with gas dynamics 
and gravity on a 64"^ N-body grid. The large wavelength 
modes of the initial density field are sampled from the ini- 
tial conditions of the large scale N-body simulations, and 
power on smaller scales is sampled from the appropriate 
CDM power s pectrum. The simulation scheme is P3MSPH 



of the form 



1 



(Evra fd 1988), the bar^ron donsit}" is a fixed fraction of the 



total Ofr^ ^r?o; and radiative rnnling is ignored 



Simulating individual clusters requires two steps: (1) an 
initial, 32'^, purely N-body simulation to identify which 
portions of the initial density field lie within the cluster 
virial region at the present epoch, and (2) a final, ef- 
fectively 64'^, three species, hydrodynamical simulation. 
In the final simulation the portion of the initial density 
field which ends up within the cluster virial region by the 
present epoch is represented using dark matter and gas 
particles of equal number, while the portions of the initial 
density field that do not end up within the cluster virial 
region by the present epoch are represented using a third, 
coUisionless, high mass, species. The high mass species is 8 
times more massive than the dark matter particles in the 
central, high resolution region. This approach allows us 
to include the tidal effects of the surrounding large scale 
structure and the gas dynamics of the cluster virial region 
with simulations that take only a few days of CPU time 
on a low end UltraSparc. The scale of the simulated re- 
gion surrounding each cluster is in the range 50-100 Mpc, 

1/3 

and varies as M^^^^, where Mhaio is approximately the 
mass enclosed within the present epoch turn around ra- 
dius. Thus, the 48 simulated clusters in our final sample 
have similar fractional mass resolution; the spatial reso- 
lution varies from 125-250 kpc. The masses of the final 
cluster sample vary by an order of magnitude. We create 
X-ray images and temperature maps for further analysis 
following procedures described in Evrard (1990). 



3. ANALYSIS 

In this section, we detail our data reduction procedures, 
placing special emphasis on the determination of statistical 
uncertainties and systematic errors in the derived cluster 
ICM gas masses. 



3.1. Modeling the X-ray Surface Brightness Profile 

We use the reduced, cluster X-ray images to constrain 
the underlying ICM distribution. Our approach is to fit 
radial surface brightnes s profiles of the cluster emission to 
the standard /3 model (Cavaliere & Fusco-Femiano 1978) 



l{R) = h 



-3/3+i 



(1) 



where Rc is the core radius and r~^^ is the asymptotic 
radial fall-off of the underlying ICM density distribution. 
Because the PSPC photon detection rate for a parcel of 
gas with constant emission measure varies onl y rn oder- 



ately with temperature above 1.5keV (see Fig. |2.l| ), the 
distribution of X-ray photons in PSPC images directly 
constrains the ICM density distribution regardless of the 
temperatu re structure (see similar discussion of the Ein- 
stein IPC: Fabricant, Lecar fc Gorenstein 198C). Fitting 



the 13 model provides a systematic approach to studying 
the ICM properties. In the case where cluster X-ray emis- 
sion is well described by (3 models — even ffatten ed 
models with axial ratios within the observed range ( Mohi 
et al. 1995) — fitting to azimuthally averaged surface 



brightness profiles recovers unbiased estimators of (3 and 
Rc as long as the point source response function (PSF) of 
the imager is properly treated. 
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Fig. 2. — In the upper panel we plot the observed surface bright- 
ness profile in A 1795 (points with error bars), the best fit single /? 
model (dashed line), and the best fit double /3 model (solid line). 
The lower panel contains residuals around the single /3 (solid) and 
the double /3 (open) fits. In clusters with central emission excesses, 
a single (3 fit produces an artificially small core radius and /3; adding 
the second /3 component results in an improved fit to the emission 
excess and the outer profile. 

Unfortunately, approximately half of nearby galaxy clus- 
ters exhibit X-ray emission which is not azimuthally 
symm etric ( |Mohr, Fabricant fc Geller 1993|, |Mohr et al' 



1995). These features include large centroid variations 
(e.g. A 754) and bimodality (e.g. A 85) which violate the 
assumptions underlying Eqn. n^. Nevertheless, we find that 
the 13 model describes azimuthally averaged cluster emis- 
sion reasonably well and produces sensible results for most 
clusters. As discussed in detail below (§3.6), we quantify 
the systematics associated with (3 model fitting to radial 
profiles by analysing an ensemble of simulated clusters that 
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exhibit X-ray morphologies similar to those of our PSPC 
cluster sample. 

The presence of cooling flows presents another potential 
problem to our choice of model. Significant central emis- 
sion excesses associated with cooling instabilities are found 
in f 8 of the clusters in our sample; they can greatly bias 
the model parameters. Clusters for which cooling flows 
present a problem to fitting are easily recognized by two 
criteria. First, the cluster must display nonrandom be- 
haviour in the fit residuals consistent with a central emis- 
sion excess. Second, the cluster must appear "relaxed", 
lacking obvious asphericity or substructure which would 

shows an example of 



3.1 



indicate a recent merger. Figure 
these criteria in Abell 1795, a symmetric cluster which dis- 
plays a significant central emission excess. The emission 
excess biases our best-fit /3 model (the dashed line) to- 
wards a small core radius and shallow profile, resulting in 
poor agreement between fit and data in the region outside 
the excess. 

In such cases we model the emission excess with a sec- 
ond (3 model and simply fit the cluster surface brightness 
profile to the sum of the two: 



(2) 



We constrain both components of this function to have 
the same /3, and we determine the distribution of the un- 
derlying ICM numerically. This method essentially decou- 
ples the inner and outer regions of the cluster, allowing 
the fit to find the transition between excess and primary 
emission on its own. An e xample of a two-component fit is 
also shown in Figure |3.l| (solid line) ; the extra degrees of 
freedom are effective in removing nonrandom trends from 
the residuals. In the process of formulating our method, 
we also considered removing the emission excess and fit- 
ting a single component (3 model to the remainder of the 
image. We found that fitting under these conditions pro- 
duces biased parameter estimates in cases where the emis- 
sion model is a perfect (3 model; specifically, the best fit 
core radius is correlated with the size of the excised re- 
gion. These systematics and associated signal to noise 
reductions led us to abandon this approach. The double 
(5 model approach is similar in spirit to this method and 
usually gives similar results, but is much more objective. 

3.2. Calculating Radial Profiles 

Producing a background subtracted cluster surface 
brightness profile appropriate for fitting requires several 
steps. First, we choose an emission center using a circular 
aperture of radius 10 pixels. We adjust the position of this 
aperture to minimize the difference between the geometric 
center and the ccntroid of the X-ray photon distribution 
which lies within the aperture. This approach converges 
even i n cases where the cluster emis sion is significantly 
skew dMohr, Fabricant fc Geller 1993| ). 

Second, we evaluate the background level in the image 
by examining annuli at large radii where the cluster surface 
brightness makes a negligible contribution. Specifically, we 
measure the background surface brightness within 5 con- 
centric, overlapping annuli with 20 pixel width and inner 
radii ranging from 140 to 180 pixels. The background value 



for each annulus is the mean of the clipped distribution of 
the ~ 2 X lO'* pixels within the annulus. The clipping algo- 
rithm excludes all pixels that are more than 2<Ti brighter 
than the median value, where ai is the uncertainty asso- 
ciated with pixel i. The clipping algorithm is iterated 10 
times, with upper and lower sigma clipping during the last 
two iterations. Tests on simulated PSPC images indicate 
that the mean value is an unbiased estimator of the image 
background. 

This method allows us to detect clusters whose emission 
contaminates the entire image, by seeing whether the back- 
ground values steadily decrease with increasing radius. It 
also guards against the possibility of localized excess emis- 
sion contaminating our result. In cases where the neigh- 
boring clusters contaminate the image background region 
in one direction (e.g. A 3558), we use half rather than full 
annuli. In the end, we measure the background in the an- 
nulus whose background value lies closest to the mean of 
the five measurements; in cases where detectable cluster 
emission essentially fills the PSPC FOV, we use the back- 
ground value in the outermost annulus. We use the scatter 
in the five background measurements as an indicator of the 
true background uncertainty. 

Lastly, we calculate the radial profile. Individual pixels 
are treated separately within 1.5 pixels of the image cen- 
troid; outside this region the profile values are averages 
over all pixels whose centers lie within an annulus of 1 
pixel width. The profile is truncated at the radius 9max 
where the signal to noise falls below a critical value; specif- 
ically, the profile is terminated when the average signal to 
noise of the last few profile points falls below '--^ 3. A few 
clusters are treated differently because of contamination 
by neighboring clusters; for example, Abell 401 is fit out 
to a radius halfway between it and Abell 399. 

3.3. Fitting f3 Models 

We find the best estimates of the parameters /q, f3, and 
Rc by minimizing the difference between the PSF con- 
volved model and the observed surface brightness p rofile 
We minimiz e using the downhill simplex method ( Press 
et al. 199^ ). In correcting for the effects of the PSF and 
finite pixel scale we employ the techniques used i n stud- 
ies of the fundamental plane of ell iptical galaxies ( Saglia 
et al. 1993| , |Mohr fc Wegner 1997| ). The PSF convolved 
mean surface brightness Ic within an annulus of radius R 
and width 2dR is 



lc{R) = 



[F{R + dR) - F{R - dR)] 



(3) 



AnRdR 

where F{R) is the integrated flux contained within a ra- 
dius R. The integrated flux is a function of the (3 model 
parameters and the PSF. 



F{R) = R / dkJi{kR)p(k)i{k) 



(4) 



where Ji is the Bessel function and p{k) and j(fc) are the 
Fourier transforms of the PSF and the one or two compo- 
nent f3 model. We use a PSF of the form 

p{k) = e-e^-^)" (5) 

where the two free parameters b and 7 describe the scale 
and shape of the PSF. We determine the best flt PSF pa- 
rameters b = 0.47 and 7 = 1.89 by fltting to artificial point 
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sources created with the PROS task rosprf at 1.0 keV (see 
Fig. |3.3| ); these parameters are appropriate for the cen- 
tral part of the PSPC field, consistent with the position of 
cluster cores in this dataset. 
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Fig. 3. — The radial nrofile (points; produced using PROS) and 
fit (solid line; see Eqn. fcl) of a 1 keV point source located in tfie 
center of the ROSAT PSFC field. The divergence of the model and 
the true PSF at radii greater than 60" does not bias estimates of (3 
model parameters. 

The Fourier transform of the one component /3 model 
I{R) is 



I{k) = 27r/oi?^ 



XJoikRcX) 



(1 + A2) 



-3/34 



where Jq is the Bessel function. This function I{k) is not 
analytic for all (3, but note that the integral is only a func- 
tion of (3 and the combination kRc- Our approach is to 
calculate I{k) over a finely sampled grid in /? and kRc and 
then interpolate on this grid when evaluating the PSF con- 
volved surface brightness I^R) (Eqn. ||). We test the ac- 
curacy of this grid by comparing the inverse transform of 
I{k) directly with I{R). The Fourier transform of the two 
component P model is simply the sum of the transforms of 
the individual components. Naturally, the PSF corrected 
surface brightness Ic{R) approaches the uncorrected pro- 
file I{R) at radii which are several times the PSF scale, 
so at large radii we use the unconvolved surface brightness 
/(i?) (Eqn. |l|) directly. 

3.4. Calculating ICM Masses 
We constrain the ICM central density using the cluster 



X-ray luminosity measured w ithin an annulus (e.g. David] 



et al. 1990, Mohr et al. 1996). Specifically, the luminosity 



within an annulus defined by a minimum radius i?_ and 
a maximum radius i?+ is 

Lx{R-,R+) ^^TT dz ne{r)nH{r)A{T)RdR (7) 

Jo J 

where ne{r) and nnir) are the electron and proton num- 
ber densities, A(T) is the radiative cooling coefficient and 



~ R^ + z^. The integral along z can be truncated 
at the cluster "boundary"; tests demonstrate that for ob- 
served values of /3 this integral is not sensitive to its upper 
limit. Using n^ir) = p{r)/ ^,f.mp and nnir) = p{r) / HH'n^p 
the ICM central density is 



Po 



' LxiR-,R+)fief^Hml{l-3f3) \ 
2nAi{Tx))R^cF{R-,R+) ) 

where F{R^, R^) is the dimensionless integral 



F(i?_,i?^ 



ds 



{l + s^ + {R+/Rc?f~'' 
[l + s^ + iR^lRcf)'-'^' 



(8) 



(9) 



and we have assumed that the ICM is isothermal with 
temperature (Tx)- For the results presented here we mea- 
sure the luminosity over the same region we fit the X-ray 
surface brightness profile (i?_ = and i?+ = dmax)- 
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Fig. 4. — The correlated errors in fit parameters /3 and Rc are ob- 
vious in this contour plot. The best fit value for A 2255 is marked 
with a star, the contours deliniate 2u and Scr confidence regions, and 
the dotted and dashed lines correspond to 90% confidence intervals 
determined by Monte Carlo runs. The uncertainties presented with 
tabulated results are derived from Monte Carlo simulations. 

In determining po we use /ie = 1.167 and [in = 1.400, 
appropr iate for a fully i onized plasma with 30% solar abun- 
dances ( |Feldman 1992 ). The ICM mass within some lim- 
iting radius is then a straightforward integral over the /3 
model density distribution using the measured ICM cen- 
tral density (Eqn. ||). 

The density distribution consistent with a two compo- 
nent (3 model is not analytic in general. We express the 
density distribution in the two component case as a ra- 
dially dependendent multiple /(r) of the broad, primary 
component: poir) = /(r)pi(r); we solve for /(r) using the 
surface brightness profiles of the two component fit Id (R) 
as a constraint: 



Id{R) 



rdr 



A((T)) 
MeMH"^p Jr Vr^~^W 



f\r)pl{r) (10) 
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Using this equation, we work our way inward solving for 
f{r) at equal intervals in r. The central density in this case 
is po = fir)piir)\^^Q. In our sample the central value of 
the multiple f{r — 0) ranges from 2.5 (Ophiuchus) to 33 
(A 2204) and averages ~10. 

3.5. Measurement Uncertainties 

A uniform, systematic analysis of the ICM properties 
in this well defined cluster sample requires careful quan- 
tificat,inn nf all known snurres nf unrprtainty We identi 



indicates recent major mergers, and, for the most 
part, they can be identified by the high xt ™ the 
fit results. We use the difference in ICM parameters 
in the two different analyses as an estimate of the 
scale of this uncertainty. The ratio of Mjcm mea- 
surements with the two centering apertures has a 
mean value of 0.997 and a standard deviation of 2.5% 
calculated over the whole sample; the largest out- 
lier is A 754, a well studied recent merger candidate 
(e.g. Fabricant et al. 198(: , Henriksen fc Markevitch 



1996|, [Roettiger, Stone fc Mushotzky 1998D, where 



multijile sources of uncertainty in the measured properties 
of the ICM, including Poisson noise, choice of cluster emis- 
sion center, PSF blurring, X-ray background subtraction, 
ICM temperature uncertainties, luminosity uncertainties, 
and (Tx) < 1-5 keV gas associated with central cooling 
instabilities. Below we discuss each in turn, and describe 
how they are incorporated into the final uncertainties for 
the P model parameters and ICM masses Micm- We also 
quantify the impact of each error contribution on the de- 
rived ICM mass, to give some indication of their relative 
importance. 

• Poisson noise: We evaluate the Poisson contribu- 
tion to the uncertainties using a Monte Carlo ap- 
proach. Specifically, we create and fit 500 indepen- 
dent realizations of each cluster image. The template 
for these artificial images is a Gaussian smoothed 
(fj = 30") version of the original cluster image. Each 
pixel in an artificial image contains a value Poisson 
distributed around the mean of the pixel value in 
the template image. When fitting j3 models to the 
artificial images, we adjust the PSF parameters to 
account for the degraded resolution in the smoothed 
template. This approach allows us to account for the 
deviations of many clusters from idealized (3 models. 
We use this ensemble of 500 fits for each cluster to 
quantify the Poisson contribution to measurement 
uncertainties. The typical Poisson contribution to 
the fractional uncertainty in Mjcm is 1.1% (sample 
median), and the largest is 3.9% in A 1689. 

The fit parameters Rc and (3 are strongly correlated. 
In Fig. ^ we present 2a and 3cr xt contours in 
the space of Rc and (3 for the single component fit 
to the cluster A 2255. Rather than present these 
plots for each cluster we note the upper and lower 
90% confidence intervals on the best fit parameters. 
Dashed lines indicate th e lo cations of these confi- 
dence boundaries in Fig. |3.4|. 



• Emission center. We use a circular aperture with 
150" radius to define the cluster emission center (see 
§3.2) before calculating the radial profile. In clus- 
ters with azimuthally symmetric emission, the emis- 
sion center is independent of the aperture radius. In 
clusters that exhibit centroid variations this is not 
the case. Calculating the profile around different 
emission centers can affect core radii, /?'s and ICM 
masses. We evaluate these effects by re-analyzing 
the cluster images using a centering aperture three 
times larger. We find that the measured properties of 
the ICM are insensitive to centering aperture radius 
except for those few clusters with grossly asymmet- 
ric X-ray images. The substructure in these clusters 



MjcM changes by 16%. 

PSF correction: Without a PSF correction, core 
radii and /3's in small core radius clusters would be 
systematically overestimated. By fitting PSF cor- 
rected (3 models to observed radial profiles, we re- 
move this bias from our P parameters and ICM 
masses. We ignore any residual errors associated 
with the minor differences between the true PSPC 



PSF and our best fit model (see Fig. 3.3 ) 



• X-ray background: Errors in the background sub- 
traction can systematically bias (3 parameters and 
X-ray luminosities. We use the variation in the mea- 
sured background within five concentric annuli as 
a measure of the true background uncertainty axB 
(see primary discussion in §3.1). We quantify the 
effects of the background uncertainty on measured 
ICM parameters by refitting the image with the 
background fixed to the preferred value perturbed 
by ±laxB and ±2axB- For each cluster we use 
the five fits to measure a maximal deviation in the 
parameters of interest (corresponding to a '--^ 4crxB 
variation in the X-ray background); we use half this 
maximal difference as an estimate of the X-ray back- 
ground induced uncertainty. In our ensemble, the 
RMS maximal deviation in M/cjv/is 3.6%, and the 
largest change is 14% in MKW 3S. 

• ICM (Tx) measurements: Typical measurement un- 
certainties in (Tx) translate into negligible uncer- 
tainties in the ICM parameters, but make important 
contributions to uncertainties in binding mass esti- 
mates. Eqn H shows that po oc y/ Lx/A where Lx is 
the X-ray luminosity and A is the cooling coefficient. 
The temperature dependence of Lx and A is approx- 
imately the same, because to determine Lx we mea- 
sure the cluster X-ray count rate in the PSPC and 
then use an ICM emission model A, galactic absorp- 
tion, cluster distance and the PSPC effective area 
tables to convert to a bolometric luminosity (PROS 
tasks). We quantify this effect by varying the tem- 
perature of the cluster A 1795 ((Tx) = 5.3 keV) from 
1 to lOkeV and calculating the resulting variations 
in Po- Above 5keV 5po ~ ST^'^^^, so a factor of two 
increase in (Tx) corresponds to a ~ 4% increase in 
Po] between IkeV and 5keV uncertainties in (Tx) 
have an even smaller effect on pQ. 

Temperature uncertainties dominate the formal un- 
certainties in the binding mass. We account for 
these uncertainties in measurements of Mjcai and 
fiCM', note that Mjcm measured within a constant 
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1/igQ Mpc aperture is unaffected by (Tx) uncertain- 
ties. 



the 90% confidence intervals. These confidence intervals 
appear with the results. 



• Luminosity measurements: We measure cluster 
aperture luminosities using the PSPC images. With 
Eqn. H and the uncertainty in the observed cluster 
count rate, we calculate the contribution to uncer- 
tainties in the ICM parameters. The typical contri- 
bution to MjcM uncertainties is 0.4% (sample me- 
dian) and the largest contribution is 1% in A 2244. 
We do not account for any errors in the PSPC effec- 
tive area tables; upcoming observations with AXAF 
will provide an excellent test of these data. 

• Core ICM with Tx < 1.5 keV: The PSPC count rate 
for a parcel of gas with constant emission measure 
is relatively insensitive to Tx for Tx > 1.5 keV (see 
Fig. 2.1); therefore, our assumption of isothermal- 



ity when we calculate the underlying density profile 
consistent with the X-ray surface brightness profile 
does not affect our results. However, ICM tempera- 
tures within the cores of clusters which exhibit strong 
emission excesses or cooling flows can be less than 
1.5 keV. In clusters where this is the case the ICM in 
the cores is essentially overrepresented with photons 
in the cluster X-ray image. Specifically, the PSPC 
count rate for a 1 keV (0.75 keV) gas clump is ~1.4 
(~ 2) times higher than the count rate for a 5 keV 
gas clump. 

In the absence of detailed ICM temperature profiles, 
we estimate the effects of these cool cores empiri- 
cally. We use the clusters with two component /3 
model fits to estimate the nature of this effect on 
MjcM measurements. Specifically, we assume that 
the secondary component of the fit is responding to 
ICM at a temperature of 0.75 keV; thus, we scale the 
best fit I2 by i to account for the factor of two higher 
PSPC count rate, and then recalculate the ICM dis- 
tribution and MicM, taking care to account for the 
change in the cluster aperture luminosity. Results 
indicate that our standard isothermal analysis over- 
estimates the true Mjcai in this case; in our sample 
of 18 clusters with double (3 model fits, AIjcm is 
overestimated by ~ 3.5% on average, with a maxi- 
mal overestimate of 14% in A 2204. The ICM radial 
distribution is less affected: (r) /rgoo is underesti- 
mated by less than 1%. 

We emphasize that the toy model presented here 
overestimates the effects of cool gas on Mjcm mea- 
surements, because we have assumed that cool gas 
extends throughout the cluster (traced by the sec- 
ondary /3 component), whereas the cool gas is likely 
confined to the core regions where the cooling time 
falls well below the Hubble time. 

We combine these sources of uncertainty as follows. 
First, these different error terms are independent of one 
another and we assume that they are each Gaussian dis- 
tributed; we combine all but the Poisson contribution in 
quadrature and call it aNP- Second, we convolve the 500 
Monte Carlo results for each cluster with a Gaussian of 
scale aNP, and integrate over the distribution to calculate 
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Fig. 5. — We plot the log of the ratiX of the inferred to the true 
Mjcm log MfQf,i) versus (Tx)for the ensemble of hydrodynamic 
simulations. Estimates of Mjcm lie systematically higher than the 
true values; the dotted line shows the median overestimate of 12% 
for the set of 144 images. 



3.6. Testing Our Technique with Simulated Clusters 

We test the accuracy of our ICM measurements using 
an ensemble of 48 hydrodynamical cluster simulations (dis- 
cussed in §2.3). The simulations allow us to examine the 
effects of (i) departures of the true ICM density distribu- 
tion from the /3 model and (ii) ICM temperature varia- 
tions within the cluster. The simulated clusters experi- 
ence present epoch mergers consistent with expectations 
for their CDM cosmogonies; X-ray morphologies of the 
simulated clusters in three of the four cosmological models 
are statistically consistent with the X-ray morphologies of 
the PSPC cluster sample, with the tCDM simulations ex- 
hibiting somewhat more substructure than observed clus- 
ters (Mohr & Evrard, in prep). We analyze images of these 
simulated clusters as though they were true cluster images 
and then compare the estimates of the ICM mass fraction 
to the real values. 

Interestingly, we overestimate the true ICM masses; the 
median overestimate for a sample of 144 images is 12%, 
and the RMS scatter around this median is ^ 9% (see Fig. 
3.5). This bias arises because the ICM distribution in sim- 
ulated clusters is not perfectly described by an isotropic, 
smooth density function. ICM density fluctuations around 
the mean value at a given radius enhance the X-ray emis- 
sivity of the gas which is proportional to the square of the 
density. We examine the ICM distribution in the simu- 
lations and measure the quantity Sicm = (p) /V {p^) ^ 1 



where p is the local ICM density; because the ICM emis- 
sivity ex scales as tx c>c p^, the overestimate in Mjcm 
scales linearly with Sjcm- We find that within rsoo the 
fluctuation amplitude is Sicm ~ 13%, in good agreement 
with the median overestimate in Mjcm- 
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These density fluctuations are presumably residual 
structure in the ICM from recent mergers; higher resolu- 
tion simulations would likely yield even higher fluctuation 
amplitudes. We are carrying out a detailed analysis of the 
ICM structure in our simulated ensemble to better charac- 
terize the source of these density fluctuations (Mathiesen, 
Evrard & Mohr in prep). In real clusters this mechanism 
and additional physics like magnetic fields and radiative 
cooling instabilities (all explicitly ignored in these simu- 
lations) could all contribute to ICM density fluctuations 
(and overestimates of Mjcm)- Quantifying the Mjcm 
measurement biases in real clusters requires detailed anal- 
ysis of high resolution, spatially resolved, X-ray spectra of 
the ICM emission or comparison of Mic m measurements 
from X-ray data and observations of the ICM's Sunyaev- 
Zeldovich effect on the cosmic microwave background. 

In summary, our ensemble of hydrodynamical cluster 
simulations indicates that (i) our analysis technique is 
likely to overestimate Mjcm, and (ii) (after accounting 



for th s bias) high signal to noise X-ray images of a coeval 
cluster population allow the determination of Mjcm with 
a typical accuracy of ~ 9%. 

4. RESULTS 

Below we present results of an analysis of our X-ray flux 
limited sample of 45 galaxy clusters. First we describe 
our analysis of the cluster X-ray surface brightness pro- 
files. Second, we present our constraints on ICM masses 
MicM and mass fractions ficM, and third, we constrain 
systematic variations in the radial distribution of the ICM. 
Finally, we compare our results to previously published 
analyses. 

4.1. X-ray Surface Brightness Profiles 

Best fit parameters for the profiles appear in Table ^ 
Columns contain the cluster name, the central surface 
brightness and core radius of the primary P model compo- 
nent, the central surface brightness and core radius of the 
secondary component (if there is one), the best fit /3, the 
reduced xt, the number of constraints in the fit, and the 
truncation radius of the fit 0„i (we also measure our aper- 
ture luminosities within 9m)- Central surface brightnesses 
are the unabsorbed values in ergs/s/cm^/D' for the 0.5- 
2.0 keV band. Core radii are in h'^Q Mpc (with go = 0.5), 
and the truncation radii 0„i are in arcminutes. Uncertain- 
ties are 90% confidence intervals with the other fit param- 
eters unconstrained. 

The X-ray surface brightness profile of each cluster ap- 
pears in Fig. ||. Each page contains a two by three ar- 
rangement of cluster panels. The clusters are arranged 
by (Tx), which facilitates the fixed radial range for each 
group of six. Cluster names appear in the upper right cor- 
ner of each panel. Surface brightness (in ergs/s/cm^/D' in 
the 0.5-2.0 keV band) is plotted versus radius in h'^Q Mpc. 
We plot measurements (points with error bars), best fit 
profile (solid line), and best fit profile without the PSF 
correction (dashed line) for each cluster. The best fit pa- 
rameters appear below the cluster name. The fit residuals 
are at the bottom of each panel. An X-ray contour plot is 
located at the lower left of each panel; contours are of a 
smoothed version of the portion of the X-ray image used 
to calculate the surface brightness profile. Contours ap- 



pear at factors of 2.5 in surface brightness, and the heavy 
contour always corresponds to 4 x 10^^'' ergs/s/cm^/D'. 

We fit two component /3 models to 18 clusters in our 
sample; these clusters exhibit trends in the residuals of 
their fits to single component P models. Some remain- 
ing clusters with single component (3 fits exhibit trends, 
but these clusters also exhibit large centroid variations or 
other evidence of recent mergers, and so we do not at- 
tempt a double component fit: (A 754, A 1367, A 2319, 
A 3266, A 3558, A 3667). Two clusters with double com- 
ponent fits also exhibit obvious evidence of recent mergers: 
A 2142 and A 85; in the case of A 85 we remove the small, 
secondary surface brightness peak to the south before cal- 
culating a radial profile. 

It is interesting to note that 29 clusters in our sample 
are not formally well fit by either a one or two compo- 
nent P model (less than 1% chance of consistency between 
model and data). This is most likely another indicator 
of the pr evalence of cluster X-ray substructure (Mohr et 
al. 1995). However, in all but a few clusters the best fit 



P model reproduces the general character of the surface 
brightness profile (no systematic trends in the residuals), 
and, as demonstrated by our analysis of the cluster simu- 
lations, this is all that is required to constrain the average 
behavior of the density profile and calculate Mjcm with 
a typical accuracy of ~9%. (Note that parameter uncer- 
tainties are determined using Monte Carlo simulations and 
not Ax^.) 

Note that the high xt fitting problem for galaxy clusters 
is similar to that of elliptical galaxies. Elliptical galaxies 
are reasonably well described by bulge plus disk models, 
but very rarely are the measured g alaxy profiles and m od- 
els consistent in a formal x^ sense (Saglia et al. 1997); in- 
terestingly, the photometric parameters derived from those 
fits exhibit the fascinating regularity relation termed the 
Fundamental Plane. 



4.2. ICM Properties 

We use the fit to the X-ray surface brightness profile of 
each cluster to constrain the ICM central density, radial 
distribution, and total mass within several limiting radii 
riirn- We use ni,„ = Ih^g Mpc and rum = ^500, the radius 
within which the mean density is 500 times the critical 
density p^rit = iH^/'&TTG. By using rsoo we are able to 
study the same portion of the virial region in each cluster; 
a drawback is that one requires a model of the potential 
to calculate rsoo, and so we also present measurements 
within a fixed metric radius to avoid this uncertainty. We 
calculate the radius rsoo and enclosed binding mass M500 
two different ways: (i) using virial scaling relations and 
(ii) using the isothermal P model. The virial scaling rela- 
tion is motivated by the simple assumptions that clusters 
have recently formed and are approximately virialized and 
self-similar; in this case rsoo oc 

and M500 oc T^/'^. 

Our simulated cluster ensemble exhibits scaling relations 
consistent wi th these assumptions, in agreemen t with pre- 
vious work ( Evrard, Metzler fc Navarro 199^ , Schindle: 
1996). On the other hand, estimates of Afsoo and rsoo with 
the P model require only that the ICM be approximately 
in hydrostatic equilibrium. Rather than using simulated 
clusters to normalize our scaling relation, we normalize 
using the P model results in the relatively relaxed cluster 
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Fig. 6h. — Ophiucus, A 1689, Triangulum Australis; Cluster X-ray surface brightness profiles: For each cluster the radially averaged, 
unabsorbed surface brightness in units of ergs/s/cm^/D'in the 0.5:2.0 keV band is plotted versus distance from the cluster emission center [in 
h'^g Mpc]. The surface brightness measurements (points with error bars), the underlying model (dashed line) and the PSF corrected model 
(solid line) are plotted. Residuals between the fit and the data (scaled by the data uncertainty) appear below each radial profile. A contour 
map of the portion of the cluster image used to calculate I{R) appears in the lower left corner; the images are smoothed to facilitate contouring. 
The emission center is marked with a star; contours appear at factors of 2.5 in surface brightness, and the heavy contour corresponds to a 
fiducial surface brightness of 4.0 X 10"^'* ergs/s/cm^/D'in each cluster. The cluster name and best fit parameters appear in the upper right 
hand corner of each plot. Both sets of parameters are listed for those clusters fit with a two component /3 model. 



A 1795. 

1 /2 

r,,o = 2.37h^,'Mpc(^^^^ (11) 

Msoo = 2.00 X lO'^h^o'Me (y^^) ^ 

If our f3 model mass for A 1795 is in error, we will still 
be studying comparable regions of each cluster; those re- 
gions would have an enclosed overdensity different than 
500 times the critical density. 

Tables || and || contain the results. The columns in 



Table || correspond to cluster name, ICM central den- 
sity po, central electron number density rig, and the ra- 
dial location (r) /r^oo of a typical ICM particle located 
within r5oo . Table ^ contains the ICM masses Mjcm and 
mass fractions ficM- The columns correspond to cluster 
name; Mjcm and ficM calculated within a limiting ra- 
dius riirn = Ih^Q Mpc; rsoo, Mjcm and //cm calculated 
using the virial relation; and rsoo, Mjcm and ficM cal- 
culated using the isothermal /3 model. All uncertainties 
correspond to 90% confidence limits (see §3.5). 

4.2.1. ICM Masses Mjcm and Mass Fractions ficM 
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Fig. contains a plot of Mjcm and ficM, calculated 
within a limiting radius 7-500, versus emission weighted 
mean ICM temperature (Tx). The left panel contains re- 
sults obtained using the isothermal /3 model to estimate 
binding masses and rsoo, and the right panel contains re- 
sults using the virial scaling relation (Eqn. |ll|). As noted 
above, measurements within rsoo sample the same por- 
tion of each cluster's virial region, whereas measurements 
within a fixed metric radius are directly comparable to 
results from previous analyses and are more nearly inde- 
pendent of assumptions about the cluster potential. 

Within rsoo the best fit relation between Mjcm and 
(Tx) is 

MicM = (1-49 ± 0.09) X 10" Mq TQ^-^>i±o-^» (12) 
where Tq is (Tx) in units of 6 kcV. The variance weighted 



RMS scatter in Mjcm around this relation is 17%; in 
calculating the variance weighted RMS, each cluster is 
weighted using its Mjcm measurement uncertainty. We 
find the best fit relation by minimizing x^, the sum of 
the squared, orthogonal deviations of each point; Vi is 
the minimum distance between point i and the fit scaled 
by one over ai, the point i uncertainty along that vector. 
We fit using only those 42 clusters with well determined 
temperature uncertainties, and we determine the 90% un- 
certainties on the fit parameters using 500 bootstrap fits. 

The 17% scatter about this relation is another indica- 
tor of the regularity of galaxy clusters; the scatter is con- 
sistent with constraints derived from analysis of scatter 
in the Size- Te mperature relation in a n X-ray flux limited 
cluster sample (Mohr & Evrard 1997) and the luminosity- 
temperature relation in a cooling-flow free cluster sample 
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Fig. 7. — We plot ICM mass Mjcm (below) and mass fraction fjCM (above) versus {Tx); these quantities are measured within the 
limiting radius rn^ = r^oo, where rsoo is calculated using the virial relation (left) and the isothermal /3 model (right; see §4.2 and Table H). 
Different symbols represent clusters with central cooling times significantly below 10'^'^ yrs (open), higher central cooling times (solid), and 
clusters with no upper limits on their (Tx) measurements (stars). 
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( Arnaud fc Evrard 1998 ). The slope of this relation is 
steep; it is statistically inconsistent with a slope of |, the 
expectation if galaxy clusters are self similar objects. We 
also examine the Mjcm-{Tx) relation using measurements 
where rsoo is determined using the virial theorem (see Eqn. 
[TT| ). In comparison to the /3 model analysis, the se dat a ex- 
hibit a somewhat smaller 14% scatter (see Fig 4.2.1 ), but 
a statistically consistent slope 1.94 ± 0.12. 

The best fit ficM-{Tx) relation within t^qq is 



fiCM 



(0.207 ±0.011) T6°-34±0-22 



(13) 



where Tg is (Tx) in units of 6 keV; the variance weighted 
scatter about the relation is 17%. These data provide ev- 
idence for a modest increase in ficM with (Tx)- The 
slope is inconsistent with zero at the 99% confidence level. 
We probe for variations in ficM by dividing the sam- 
ple at (Tx) = 5 keV; the weighted mean of the 17 low 
(Tx) (28 high (Tx)) clusters is 0.167±0.008 (0.212±0.006), 
corresponding to a statistically significant difference of 
0.045±0.010 (uncertainties are 90% confidence intervals 
determined using RMS scatter of measurements about the 
mean). The mean ficM of the low (Tx) sample rises to 
0.176±0.006 if A 1060 is excluded, bringing the difference 
in the mean of the two samples to 0.036±0.001. An anal- 
ysis of ficM where r5oo and M500 are calculated using a 
virial scaling relation suggests a slightly more rapid in- 
crease of ficM with (Tx). The best fit slope is 0.41±0.16, 
and ficAi is larger in hot clusters than in cool ones by 
0.052±0.009. We caution that the virial scaling rela- 
tion assumes self similarity in the cluster population, but 
the measurements presented above demonstrate the ICM 
properties are inconsistent with self similarity. Therefore, 
it is a distinct possibility that the virial scaling relation 
masses M500 are systematically in error. 

We have examined the same relations within a fixed ra- 
dius of Ih'^Q Mpc. The best fit relations between Mjcm, 
fiCM and (Tx) are 



MicM = (7.45 ± 0.49) x 10^^ Mq Tgi-23±o.i7 ^^^^ 
fiCAi = (0.184 ±0.009) 



(0.184 ± 0.009) Tg°-i^±" " (15) 



where Tg is (Tx) in units of 6 keV. As expected, the 
MicM-{Tx) is shallower than the corresponding rsoo re- 
lation presented above. The fiCM-{Tx) relation is flatter 
and consistent with no increase of ficM with (Tx)', this 
is primarily caused by a reduction in ficM in the higher 
(Tx) portion of the sample, behavior that implies an in- 
creasing ficM with radius. 




4 



6 8 10 



Fig. 8. — We plot ICM mass Mjcm (below) and mass fraction 
fiCM (above) versus (Tx); these quantities are measured within 



the limiting radius rii„ 



IfegQ Mpc. Clusters with different cen- 



al cooling times are represented with different symbols as in Fig. 



In addition, the data suggest a different Micm-{Tx) 
zeropoint in short and long cooling time clusters; the ze- 
ropoint for the 20 short cooling time clusters is 8.19^°'^^ x 
lO^'^ Mq and for the the 22 long cooling time clusters is 



±0.43 

^±0.92 



X 1013 M 



This corresponds to AM, 



ICM 



6.29 

1.90±° '32 X 1013 Mq, a 3.3cr or -25% offset. An offset 
of similar scale appears in measurements made within a 
0.75h^o Mpc aperture; AMicm = 1.32±"-65 x lO^^ M© 
which is a 3.2a or ~25% offset. At the larger radius r5oo 
the offset is AMicm = 2.73^^-^'^ x lO^^ Mq which is a 
2.3cr or 18% offset. 

This correlation between Mjcm and central cooling 
time could be because, at a particular mass scale or (Tx), 
clusters with higher gas content are more likely to expe- 
rience central cooling instabilities; in that case we need 
to understand what mechanism drives the differences in 
gas content. Our numerical cluster simulations which in- 
clude gravity and gas dynamics provide no evidence for 
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large Mjcm offsets like those observed. Alternatively, one 
might suggest that we are seeing the results of the ICM 
"flowing" into the cluster core; the scale of the mass off- 
set requires an average flow rate of ^--^2000 Mo/yr over 
the age of the universe, which makes this explanation un- 
likely. On the other hand, the higher Mjcm in clusters 
with shorter central cooling times may be an artifact; cen- 
tral cooling instabilities may enable the formation of a 
multiphase medium and significant ICM density and tem- 
perature variations rather than the smooth density distri- 
bution and isothermal ICM implicit in the /3 model anal- 
ysis. These density variations enhance the X-ray emissiv- 
ity, causing an increase in the X-ray luminosity and a bias 
in our estimates of the central density po (see Eqn. ||). 
Quantitatively the offsets in Mjcm could be explained by 
a multiphase medium in short central cooling time clusters 
if that medium enhanced the X-ray luminosity by ~ 50% 



the core and asymptotic behavior of the profile to be fit 
simultaneously; the gross trends in the residuals disappear 
and the goodness of fit parameter falls by a factor of 
5. As discussed in §3.1, adding the second component is 
similar in spirit to excluding the core region from the fit; 
however, results of fits to core excised clusters depend on 
the scale of the region excised, exposing the fit parameters 
to the subjective whims of the observer. We choose double 
P model fitting because it is more objective, and because 
our data demonstrate that the extra component allows the 
core and the asymptotic behavior of the profiles to both 
be fit simultaneously. 

We also probe for trends in the ICM distribution with 
(Tx) using Peffi a nonparamctric measure of the slope of 
the surface brightness profile. Specifically, 



on average. 



4.2.2. ICM Radial Distrihutic 



log(/./4) 1 
\og{Ao/A{) ^ 2 



(17) 



We also probe for systematic differences in the radial dis- 
tribution of the ICM in high and low {Tx) clusters. The 
presence of differences could be an indication that galac- 
tic winds have contributed significant energy to the ICM 
(IWhitc 199l| , iMetzlcr fc Evrard 1994 [Metzler fc Evrard 
1998 ). In past studies, discussion has focused on the value 



-3/3, the asymptotic slope of the ICM; we take a new 
approach which incorporates the entire ICM distribution 
within the radius rsoo. Specifically, we calculate the aver- 
age radial location (r) / rsoo of an ICM particle within the 
same physical region of each cluster, where 



operationally, we choose an area Aq, find the isophote 

10 which encloses that area, and then find the isophote 

11 which encloses the area Ai = Q.IAq. This approach 
yields the effective ICM falloff —Weff over a particular 
annulus determined by Aq. Fig. (upper panel) con- 
tains a plot of f3eff measured at rsooo versus (Tx) for the 
same sample; the data provide no indication that lower 
(Tx) clusters have a more extended ICM. Because Peff 
changes with radius in a cluster (small in the core and 
steepening with radius) it is important that it be mea- 
sured in the same physical region in each cluster (rathe r 
than at a particular isophote; see Mohr fc Evrard 1997 ). 
The results displayed here are calculated at a radius r = 

Differences in the ICM radial distribution with (Tx) will 0.75h^^{{Tx) /lO keV) Mpc, corresponding to rgooo and 
appear as trends in (r) /r^oo versus (Tx)- Ao 

Fig. (lower panel) contains a plot of (r) /rsoo versus each cluster image. 
{Tx) for our cluster sample. Symbols code central cool- 
ing times as in other figures. The distribution of (r) /rsoo 
in this flux limited cluster sample provides no indication 
that the ICM is more extended in low mass clusters. These 
data stand in contrast to previous analy ses of large Ein- 



/q""" (Pr rpjr) 



(16) 



a radius at which P^ff can be measured in 



4.3. Comparison to Previous Results 



stein IPC samples (e.g. Mohr et al. 1995) which indicated 
a tendency for low (Tx) clusters to have lower /3's (more 
extended ICM distributions). 

We have analyzed the entire sample with a single /3 
model fit, and we recover a trend for lower (Tx) clus- 
ter to have a more extended ICM. Therefore, the root of 
the discrepancy between our results and previous results 
lies primaril y in the treatment of the central emission ex- 



cesses. Fig. 3.1 demonstrates the differences between sin- 



gle and double component /3 model fits to the X-ray data 
on A 1795, a cluster with a central emission excess. We 
plot the best fit single component model (dashed hue), the 
double component model (solid line) , the residuals around 
the fit (single: solid points, double: hollow points), and we 
list the best fit parameters and reduced x^- The behavior 
in A 1795 is generic to clusters with central emission ex- 
cesses; the best fit single /? model has a small core radius 
and correspondingly low (3 that grossly misrepresents the 
asymptotic behavior of the surface brightness profile. The 
two extra free parameters in the double /3 model fit allow 



We compare our Mjcm and fjc m measurements to pre- 
viously published results. Rather than present detailed 
notes on each cluster, we focus on two trends apparent 
from the comparisons. First, our measurements in the 
clusters fit with single (3 models are generally consistent 
with previously published results. Second, our measure- 
ments on clusters with central emission excesses tend to 
produce lower ICM masses Mjcm than previously pub- 
lished results. In individual cases of disagreement in /3 
parameters, we invite the reader to examine the X-ray 
surface brightness profiles and contour plots (Fig, n j as 



imele, Ncsci fc Trevesc 1997], David, Jones fc Formar 



a way of gauging the plausibility of our fit ( Allen et al 
1996 L iDardelli et al. 1996L i Brecn et al. 1994 |Briel ct al 
19921 IBriel fc Henry 1996|, iBuote fc Canizares 19961, ICn" 



19951, IDeU'Antonio, Geller fc Fabricant 19951 Fabricant et 



al 19891, IHenry et al. 1993|, penriksen fc White 1996|, pen 



riksen fc Jones 1996, IMakTno 1994, [Neumann fc Boehringc: 



fcFabian 1995) 



19951, Pislar et al. 1997| , [Sarazin fc McNamara 1997] , [White 
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Fig. 9. — We plot the typical location of a gas particle within rsoo 
(below) and /9e//i ^ nonparametric estimate of the radial falloflf of 
the ICM calculated at ~ rsooo (above). As before, symbols denote 
different central cooling times. Neither dataset provides evidence 
for a more extended ICM in low (Tx) clusters (see discussion in 
§4.2.2). 

We have examined the issue of emission excesses in 
great detail; we analyze all clusters with central emission 
excesses using a single and double component /3 model. 
In the median, the double (3 model analysis provides an 
MicM which is 7% lower than the single (3 model analysis 
(maximum reduction is 20% in A 3526). This results be- 
cause the emission excess biases the parameters Rc and P 
low in the single component fits. The works which re- 
ported higher ICM masses generally used a single-beta 
model approach and ignored the central emission excesses; 
in these cases we were able to demonstrate consistency by 
showing that we could match the author's result using a 
single-beta model ourselves. A few authors took the ap- 
proach of excising the central region before fitting the clus- 
ter profile; these works generally show better agreement 
with our double-beta model results. 

We also found four papers which do not follow the trends 



outlined above. Three of these ('Durret et al 1994 , Myers 
et al 19971, iFukumoto fc Ikeuchi 19921 ) employ" very dif- 



ferent analyses, making it difficult to discern the source 
of the disagreement. The recent dcprojcction analysis of 
207 c ooling flows using IPC data (White, Jones & Forman 



1997) presents ICM masses and mass fractions for many 



of our clusters and seems to employ a comparable method, 
but we found poor agreement between our results. There 
are 18 clusters which we both measure out to lh~^ Mpc; 
out of this set, there are 9 clusters for which our uncer- 
tainty ranges in Micm do not overlap. Strangely, there is 
no evidence of a systematic difference in our results; when 
plotted against one another the scatter is just much larger 
than the error bars. Analysis of the ICM mass fractions 
fiCM is even more discouraging, with only two clusters in 
agreement. Here, however, there is a systematic difference 
in that the fractions of WJF are almost all much lower than 
our own. The worst offenders are clusters for which WJF 
uses older X-ray temperature determinations (Cygnus A, 
A 2029, A 2142, A 2319, A 401 and A 496) or which have 
no measured velocity dispersions. The source of this dis- 
crepancy is still unclear, but seems to be due at least in 
part to the large uncertainties in the galaxy velocity dis- 
persions agai , which sets the depth of the potential well in 
their deprojection technique. As already mentioned in the 
introduction, calculating ICM masses and mass fractions 
docs not require assumptions about the detailed shape of 
the cluster potential, and so we make none in our analysis. 

To conclude, we reiterate that because the ICM distri- 
bution follows directly from the X-ray surface brightness 
profile, a model which provides a better fit to the data will 
yield a more accurate estimate of Mjcm (see §3. Exten- 
sive testing of the double beta model approach convinces 
us that the technique produces a significantly better de- 
scription of the data, allowing simultaneous fitting of the 
core emission excess and the asymptotic behavior of the 
gas at the limiting radius of the data. 

5. DISCUSSION 

We analyze the X-ray and ICM properties of an essen- 
tially X-ray flux limited sample of 45 galaxy clusters. The 
fundamental observations are archival ROSAT PSPC im- 
ages and published {Tx) measurements; we describe our 
techniques in §3 to facilitate reproducibility. Combining 
our data and an ensemble of numerical cluster simulations, 
we quantify the uncertainties in the measured quantities. 

In addition to constraining the properties of the ICM in 
nearby clusters (see discussion below and §4), we include a 
set (see Fig. ^ of cluster X-ray surface brightness profiles 
in physical units along with contour plots of the cluster 
X-ray emission to help the reader evaluate the goodness of 
fit, the importance of substructure and the overall plausi- 
bility of the analysis in each cluster. While a few clusters 
in our sample are undergoing major mergers, making our 
analysis suspect, it should be remembered that our anal- 
ysis of an ensemble of simulated clusters indicates that 
MjcM can be estimated with a typical accuracy of ^--^9% 
in coeval populations with similar morphologies to those 
observed (i.e. even with merging clusters included). 

5.1. ICM Regularity 

The observed scatter around the Mjcm-{Tx) and ficM- 
{Tx) relation is 17%; correcting the scatter for the 9% 
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observational contribution (derived from numerical simu- 
lations) reduces the scatter to 14%. This scatter is consis- 
tent with predictions made using the scatter around th e 
cluster Size- Temperature relation ( Mohr fc Evrard 1997 ). 
The properties of this X-ray flux limited cluster sample 
provide further evidence tha t galaxy clusters are, on aver- 
age, r egular obj ects (see also Mohr fc Evrard 1997 , Arnaucj 
& E vard 1998 ). Interestingly, "exceptional" clusters do 
exist; A 1060 has an Mjcm and ficM whi ch are low rel- 
ative to clusters with similar (Tx) (see also Lowenstein & 
Mush|3tzky 1996]). 



5.2. Self- similarity and Evidence for Galaxy Feedback 

Measurements of Mjcm at r^Qo provide evidence that 
the ICM content of galaxy clusters is not self-similar. 
Specifically, the slope of the Mjcm-{Tx) relation at rsoo 
is 1.98±0.18 (90% confidence unless otherwise noted), dif- 
ferent from the self-similarity expectation of | at 4.3cr. 
We also probe for trends in ficM with (Tx)- The best fit 
fiCM-{Tx) relation has a slope of 0.34±0.22, and is incon- 
sistent with no variation at 2.5cr. If the steep slope of the 
MicM-{Tx) relation were entirely caused by an increase 
of fjcM with (Tx), the required ficM-{Tx) relation slope 
would be 0.48, well within the 90% confidence interval of 
our measurements; therefore, it is possible that the steep 
MicM-{Tx) slope is due solely to an increase of ficM with 
(Tx). We also examine trends in ficM by comparing clus- 
ters with {Tx) above and below 5 keV. We find that the 
28 clusters with (Tx) > 5 keV have //cm=0.212±Q.006, 
0.045±0.010 higher than the mean of the 17 low (Tx) clus- 
ters. Thus, these data indicate that ficAi increases with 
(Tx), and that this increase may be responsible for the 
steepness of the Mjcm-{Tx) relation. 

This departure from self similarity in the ICM properties 
of clusters could be a reflection of a lack of self similarity 
in the underlying dark matter distribution; the "universal" 
NEW density proflle predicts significant departures from 
self similarity in the cluster population driven by earlier 
formation, on average, of low mass clusters when comp ared 
to high mass clusters ( [Navarro, Frenk, fc White 1997| ). On 
the other hand, the ICM is sensitive to a wider range of 
physics than the dark matter, and this sensitivity might, in 
part, be responsible for the observed departures from self 
similarity. Energy injection into the intergalactic medium 
by supernovae driven galactic winds is one possible mech- 
anism; a trend of increasing ficM with (Tx) is consistent 



with some models of galaxy feedback ( Metzler & Evrard 
19981, [Tozzi, Cavaliere fc Menci 1998|). Another possibiL 



ity is that th e efficiency of galaxy formation varie s with 
cluster mass ( David ct al. 19901 , Arnaud et al 1992). 

Interestingly, our extensive dataset provides no evidence 
that the ICM is more extended in low (Tx) clusters. 
Specifically, there are no apparent trends in the typical 
radial location (r) /rsoo of an ICM particle within ^500, a 
non-parametric estimate of the radial fall-off of the ICM 
density /?e// at rsooo, or /3 measured fitting to the whole 
surface brightness profile. Previous analyses of large sam- 
ples have relied pr ima rily on single /3 model fitting, which 
we show (see Fig. 3.1 and §4.2) returns highly biased esti- 
mates of P and Tc in clusters with large, central emission 
excesses. This bias, combined with the greater prevalence 
of low (Tx) clusters with central emission excesses, leads 
to an artificial trend in /? or (r) /r^oo with (Tx). We re- 



cover this trend when we analyze our entire sample using 
solely single f3 model fits. Models of galaxy feedback gen- 
erally produce a more ex tended ICM in low mass clusters 
( Mctzler fc Evrard 1998| ). While our data do not rule out 
a weak trend in ICM extent, they provide no suggestion 
that one exists above 2.4 keV. 

5.3. Constraints on f2jv/ 

Under the fair sample hypothesis, baryon fraction mea- 
surements within cluster virial regions can be used to de- 
rive an unbiased estimate of the universal baryon fraction. 
By combining measurements of ficM in clusters (a lower 
limit on the cluster baryon fraction) with nucleosynthe- 
sis constraints on the baryon to photon ratio, we place 
an upper limit on the cos mological density p arameter for 
clustered matter V,m (e.g. White et al. 1993). 

Because more massive clusters are less affected by pro- 
cesses like galaxy feedback which might enhance ICM de- 
pletion within the cluster virial region, we use the mean 
fiCM calculated in the (Tx) > 5 keV subsample of 27 
clusters: //cAf =0.212/i5q'^^^ (see §4.2.1). We apply two 
corrections to this mean ficM'- (1) a correction for ICM 
depletion within rsoo (driven by shocks during cluster for- 
mation), and (2) a correction for ICM clumping which 
causes an overestimate of ficM- We also use the numer- 
ical simulations to estimate ICM depletion within rsoo; 
our ensemble of 48 simulations indicates that, on aver- 
age, ficM within rsoo is 12% lower than the cosmic value. 
As discussed in §3.5 and §4.2.1, ICM density variations 
or clumping enhance the X-ray luminosity and bias ficM 
estimates high. Analysis of our numerical simulations in- 
dicates that, on average, estimates of Mjcm and ficM are 
too high by 12%. There is no guarantee that the struc- 
ture responsible for this overestimate in our simulations 
mimics the structure present in real clusters. Hot galactic 
winds, magnetic ffelds, and radiative cooling instabilities 
could well enhance density structure; in fact, larger varia- 
tions in the ICM density distribution, and therefore larger 
systematic errors in Mjcm, ar e possible if the ICM h as a 
strongly multi-phase structure (Gunn fc Thomas 1996, Na- 
gai, Sulkanen & Evrard in prep) . Comparisons of cluster 
X-ray emission and the Sunyaev-Zel'dovich effect decre- 
ment in the cosmic microwave background should further 
constrain the density and temperature variations in the 
ICM. 

Correcting the mean ficM for clumping and depletion 
(these corrections cancel), we get our best estimate of the 



cosmic ICM mass fraction: ficM = (0.212 ± Om&)h^ 



'50 



where the uncertainty represents the 90% confidence range 
from observational uncertainties (§3.5). Constraints on the 
clustered mass density parameter £7 m are then 



ICM, 



- (0.36±0.01)e,y,/i5o'^' (18) 



where we h ave used (0.077 ± 0.001)/i5(; ( [Buries fc 



Tytler 1998 ). The ~ 3% uncertainty reflects the observa- 
tional uncertainty of the mean ficM] Qsys is a scale factor 
that represents systematic effects which we haven't ad- 
dressed. Our analysis of an ensemble of 48 numerical clus- 
ter simulations suggests that after correcting for depletion 
and gas clumping, Qgys = 1- Deviations of Qsys from 1 
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would reflect systematics which appear in real clusters but 
not i n simulated clusters. Some possibilities are (1) very 



mass analyses of even larger cluster samples a nd observa- 
tions from the well calibrated AXAF satellite ( Eisner et a] 



steep temperature profiles at rspp and (2) a very strong 1994 ) will further constrain Q 



multiphase medium extending throughout the cluster. Re- 
cent comparisons of lensing mass estimates and isothermal 
P model binding mass estimates also suggest Osys ^ 1; 
the mean ratio of isothermal /? model to lensing b inding 



Jsys- 



mass es is 1.04±0.07 for the CNOC cluster sample (Lewi; 
et al ] 998 ) . Additionally, cluster ICM masses Mjcm have 



an associated systematic uncertainty due to uncertainties 
in the absolute flux calibration of the ROSAT PSPC; we 
estimate that the PSPC absolute flux calibration is good 
to ^ 15%, implying that ICM masses are uncertain by ap- 
proximately 7.5%. Combining these estimates of system- 



atic uncertainties yields O 



sys 



1 ± 0.10. Future binding 
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Table 2 

Parameters for X-ray Surface Brightness Profiles 
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ICM Density and Distribution 
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C C 1 u.uuy 
■'■0.010 


A3112 


'^•^lato^' 


■2 


7.040:t|e- 


•26 


0, 


r770.015 
•3' '0.015 


A 3158 




-3 


8.08S]:ifc- 


■27 


0, 


c: QqO.007 
••J ''•Jo. 007 


A 3266 


OT 

^•'j8o;o7'^" 


-3 


4.63S]:l|c- 


■27 


0, 


gQQO.006 
•"^*J0.006 


A 3391 


9 7')0.30 


-3 


r, 01 0.59p 


■27 


0, 


fil 40.010 

•o^i^o.oio 


A 3526 


6.8lS]:I«c- 


-2 


1 oo0.15_ 


■25 


0, 


trqi 0.010 

■•jyio.olo 




9 4q0.34 
^•^^0.27'- 


.Q 


4 Sfi0-66p_ 
^-oOo.SS*- 


■ 97 


0, 


finso.oio 

•"U»o.oi3 


A 3558 


4.95g:iie- 


-3 


9.65S]:i|c- 


■27 


0, 


f;nQ0.006 
•OU»o.o06 


A 3562 


5.880;ie- 


-3 


1 1 c;0.08 
-'^•^•JO.05'- 


26 


0, 


6260 002 
•0^00 002 


A 3571 


7 090. 23p 
' -^^0.29*^ 


3 




26 


0, 


C770.007 

■3' '0.007 


A 3667 




-3 




■27 


0, 


pif;0.003 
•"^"0.003 


A 4038 


l-43g:Jie- 


-2 


2.78°:ie- 


26 


0, 


c;q90.013 
■'3^^0.013 


A 4059 


1 iq0.09p. 

-'-•-'-^0.06"= 


-2 


9 Ql0.18„_ 

^•^-•^o.is*^ 


■26 


0, 


5860-007 

•3°'J0.006 


0745-19 


6.840:ie- 


-2 




■25 


0, 


c:cq0.004 
•'J"'J0.004 


AWM 7 




-2 




■26 


0, 


c:(;f;0.006 
•'J""0.006 


Cygnus A 




-1 


9 C90.71„ 
^•"J^0.34<- 


■25 


0, 


fil 40-021 
•01^0.021 


MKW3S 




-2 


"^•"^^0.20*^ 


26 


0, 


cyqO.OlS 
■^'^0.015 


Ophiuchus 




-2 


3.08;]:ic- 


26 


0, 


r^RO.OOQ 
•"J^Oo 009 


Tria jVusI 




■3 




■2(i 


0, 


■•""(1.(114 



26 



Table 4 

ICM Masses and Mass Fractions 



Cluster 



MiCM 



Ih, 



-5/2 



5o' Mpc 

fiCM 



"50 



-3/2 



Ko Mpc 



Virial Relation at rsoo 

MjCM fiCM 



-5/2 



"50 



-3/2 



"50 



rii 



Mpc 



(3 Model at r5oo 

MjCM fiCM 



A 85 
A 119 
A 262 
A 401 
A 426 
A 478 
A 496 
A 644 
A 754 
A 780 
A 1060 
A 1367 
A 1651 
A 1656 
A 1689 
A 1795 
A 2029 
A 2052 
A 2063 
A 2142 
A 2199 
A 2204 
A 2244 
A 2255 
A 2256 
A 2319 
A 2597 
A3112 
A3158 
A 3266 
A 3391 
A 3526 
A 3532 
A 3558 
A 3562 
A 3571 
A 3667 
A 4038 
A 4059 
0745-19 
AWM7 
Cygnus A 
MKW3S 
Ophiuchus 
Tria Aust 



7.95°; 
5.42° 
2.44g 
1.02° 
9.33[] 
1.020 
4.93g 
7.85° 
8.72[; 
5.78° 
1.99° 
3.71° 
7.588 
8.31U 
1.181] 
7.968 
1.068 
3.638 
3.998 
1.298 
4.778 
1.178 
7.498 
7.138 
9.178 
1.108 
4.998 
5.848 
6.898 
8.408 
4.838' 
2.668' 
4.578 
7.378 
4.758 
7.938 
8.168 
3.488 
4.248 
1.208 
3.848 
7.708 
3.818 
9.858 
1.088; 



.866+13 
.866+13 

■8ie+14 
:8?c+13 
.8ie+14 
.866+13 
■i8e+13 
■1^6+13 
.85C+I3 
.856+13 

■866+13 

■J86+13 

08 
.08 

"2c+14 



c+13 



.03 

:88c+i3 

•8}c+14 

580+13 
:8ic+i3 
:88e+i4 

■836+13 

8^e+14 

i6+i3 

■o|6+13 
•}8c+13 
8ic+14 

\te+n 
■156+13 

.856+13 
■\lc+13 
•}8c+13 
l\c+13 
896+13 
:i8e+13 
•8ie+13 
876+13 
.836+14 
8|c+13 
■|ie+13 
■186+13 

12c+13 

8ie+14 



0.1948:88? 

0.1588:818 
0.1688:8}?, 
u-J^y^Q 008 

n r,r)70.022 
^•^^"0.021 

U.ZU^o 006 

0.1848:888 
0.1698:888 

n 1 790.007 
U--"-' ^0.007 

n 1 c«o.oo8 

'J--'-°Do.008 
PI r)oc0.008 
U.U»Oq 007 

0.1798:881 
0.1828:818 

n 1 tcnO.005 
U.iOUo.oo5 
n 1««0.070 
U.iDDo.o69 
A 1 070.005 
U.J^y '0.005 
n 17CO.OO5 
'-'•J- '^0.005 

U.J-^Oq 019 
n 1970.021 
".J^^' 0.020 
A 991 0.007 
u.zzig 006 

1 c;i 0.010 

'J-J-3J-0.009 

".^J^J^o.oii 

•J. -"-0^0.047 
A 1 r aO.030 
".^^^'-'0.030 
A iee0.003 
U.iDDo.o03 
n 01 q0.004 
U.ZiOo 004 
A 9r)f;0.015 
U./:UDo.oi5 
A 9A9O.OI5 
U.ZUZq 

A iocO.014 
U.iSOo.oi4 

A nncO.OlG 
U./UDq 016 
n 1 4Q0014 
".^^^0.013 
A 199O.OIO 
^.^^^0.009 
A 1700013 
".J^ ' ^^0.015 
A 900O.OO6 

U.ZZOg 006 
A 944O.O34 

".^^^0.034 
A i7f;0.006 
"J -J- "Jo. 006 

0.2098:811 
0.1798:811 

A 1700010 

U. ZZUq 027 

0.1378:881 
o.2298:81i 

".J^ 'O0.016 

A IQ4O.OIO 

^•-'■<5%.010 

A I70O.OI5 

'J-^' "^0.013 



1.873 
1.826 
1.177 
2.185 
1.908 
1.983 
1.500 
1.947 
2.211 
1.478 
1.335 
1.419 
1.903 
2.173 
2.410 
1.752 
2.237 
1.335 
1.536 
2.234 
1.609 
2.275 
2.021 
2.049 
2.078 
2.290 
1.439 
1.644 
1.778 
1.888 
1.810 
1.427 
1.591 
1.810 
1.478 
1.992 
2.006 
1.378 
1.536 
2.211 
1.498 
1.933 
1.419 
2.374 
2.337 



1.838:8ie+14 



1.418: 



3.088 
3.068 
2.088' 
2.248 
8.328' 

1.778: 
2.918: 

8.798 
2.768 
6.508 
1.768 
2.398 
3.338 

1.488: 

2.788' 

5.078: 

6.768 
3.708: 



A IQ9O.OO5 

•J.J^y^o.oos 

U.lDUo.oiO 
A 1 oaO.005 
U.lOUo 005 

A 9A9O.OO9 
U.ZUZo 009 

c+14 0.2068:8o9 



03'= 

856+14 

llc+13 
8ic+14 

08 
07' 

856+14 

}?e+13 

88c+i4 

31 



31 



?.408 

.528: 
,898: 
.028: 
328: 
958: 



e+14 
lie+13 

lie+13 
lic+13 
■8ic+14 
llc+14 
??c+14 
8ic+14 

06^ 

oe*- 

:^}c+i3 

i5c+13 
S^6+14 

18e+13 

22, 



ic+U 



-14 



7.738' 

1.148: 



1 4Q9-' 
2.11 



23"- 

ic+14 
18C+14 
856+14 
llc+14 
ic+13 
8?e+14 
06e+14 



0.06' 
0.09^ 



1.228: 

4.388 

9.138: 
1.818: 

9.048 
1.9981 
2.468: 
5.488: 
7.688 

3.318: 

6.298 

2.198: 
6.138: 

2.758' 



1+14 
8|e+14 
iJe+13 
ie+13 

03 c 
03*= 
65r 
65'= 

856+14 

8?e+14 

:^«c+13 



286+13 

ie+14 

ic+13 
llc+14 
|?c+13 
i5c+14 



A IA7O.OO6 

U.J-y' 0.006 

A 1 7AO.OO4 

U.J- 'Uo.004 

A icco.ooe 
U.1DO0.006 

A ioc;0.020 
U.IOO0.02O 

0.1878:881 
0.0808:881 

n 1 C7O.OO6 

'J.-LO'o.ooe 

n -|7f!0.009 
L'.J^ '"0.009 
A icaO.008 
U.iDUo.008 

0.1648:888 
0.1908:881 

n 171 0.005 
U.J- '-"-0.005 
A I4f:0.013 
^.-'-'100.012 
A 19QO.O2O 

u. Jzyo.020 

A 99Q0.004 
U.ZZ»o.o04 

n 1 '?qooo5 
u.ioyo.oo5 

A onfiO.ois 

U.ZUOo.014 
A -1 coO.046 
U.lOOo.045 

^.-^"^0.032 
A 1 70O.OO5 

u.-L '80.005 

A 99f;0.007 
U.ZZDo.007 

0.1788:811 
0.1778:811 
0.1838:811 

A 9-|f;0.017 
U.^-LOo.017 
A 141 0.012 
U.J-^-Lo.012 
A 1 A4O.OO5 

A -1 cf;0.006 
U.-LOO0.OO7 
0.005 



:+i4 0.2108:88^ 

:+13 



3.378:lie+14 



A -|q')0.026 
U.iaOo 026 

A 1 70O.OO7 

U.-L ' -50.007 

A onqo oii 
u.zu»o.oii 

A I44O.OO9 

A i4f;0.008 
".-'-'±00.008 
n 91 10.026 
'-'.Zllo.026 

0.1298:888 
0.2098:811 
0.1488:811 
0.1428:811 
0.1828:811 



1.719 
1.635 
0.998 
1.934 
1.849 
1.901 
1.369 
1.805 
1.953 
1.464 
1.273 
1.216 
1.706 
2.056 
2.221 
1.752 
2.128 
1.276 
1.458 
2.206 
1.493 
1.979 
1.782 
2.008 
2.113 
1.911 
1.290 
1.412 
1.632 
1.801 
1.511 
1.227 
1.378 
1.524 
1.157 
1.776 
1.673 
1.156 
1.313 
1.940 
1.400 
1.523 
1.219 
2.270 
2.340 



1.648: 
1.208: 

2.428' 
2.608' 

2.018-846+14 



8^e+14 
856+14 
816+13 
8?6+14 



2.148: 

7.428' 

1.628: 

2.458: 
8.708' 
2.628' 
5.118; 

1.538: 

2.238 
3.048' 

1.488: 

2.638' 

4.818: 

6.368' 
3.658' 

7.718: 

2.938' 
1.618: 
1.978: 
2.368' 
3.018: 
6.788' 
9.308 

1.348: 

1.998: 

9.248' 
3.558' 
7.428' 

1.418: 

6.058' 
1.728: 

1.868: 

4.278' 

6.198: 
2.818: 
5.818: 
1.508: 

4.988' 
2.628' 
3.388' 



1.04"- 

8te+14 
l?e+13 

03, 
03"- 
44 



-14 



1^6+13 

8ic+13 
ic+13 

!]|6+14 

:8ic+i4 

8|c+14 

836+14 

8;5c+14 
?86+13 
io6+13 
ic+14 
l|e+13 
lle+14 
ic+14 
180+14 
836+14 
I06+I4 
l|o+13 
ie+13 
8ie+14 

09, 
09"- 

i86+i3 

8|6+13 

:l6+i3 
856+14 

^|e+13 
836+14 
8?e+14 
8|e+13 
l?e+13 
116+14 
lic+13 
l|o+14 
116+13 
186+14 

1^6+14 



A O0'?0.009 
U.ZZOq 009 

0.1908:811 
0.1688:818 
0.2478:811 
0.2198:811 

A 9I4O.OI2 

"•^±%.011 

0.1998:881 
0.1908:811 

0.136 



e+14 0.2268:i^8 



0.1918:88^ 
0.0878:888 

n iqro.009 

U.J^yOo 009 

A 919O.OI2 

".^-■-^0.012 

A 1 77O.OI9 

U.J-' '0.019 

A 1 01 0.083 
U.iyio.080 

0.1908:881 
0.1888:818 

A iRoo o^e 

U.IDUq 025 
A 141 0.025 
U.±'*±0.024 
A 997O.024 
U. ^^'0.017 
A 1 cqO.006 

u.ioyo.oo6 
0.2608:81 

0.1968:88?, 
0.1688:81 

A -I79O.OI6 

'-'.-'-'^0.016 

A 9070.027 

U.zy/g 027 

A 917O.O2O 

".^± '0.020 

A 997O.O27 

U.^^' 0.027 

A 91 9O.O2O 

^.^-•-^0.019 

0.021 

0.1848:1 

A 1 09O.OI2 
U.iOZg 012 
A 1 QCO.OIS 

U.iyOg 022 

o.2738:81i 
0.2698:81^ 

A 91 10.013 
'-'.^-L±0.013 

o.2748:81i 
0.1908:818 
0.1888:81^ 

A 9G4O.O34 

U.ZD% 033 
A 14«0.010 
U.J*Do.010 
A 9090.032 
U.zyZg 032 

0.1898:811 

A 1 C4O.OI6 

U. -1-0%. 016 
0.181o:o22 



-14 0.2348 



